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Abstract 



We consider various special cases of gluon condensates at finite temperature. The 
gluon condensate for an ideal gas of gluons with a given vacuum expectation value is 
introduced for the sake of comparison with that calculated using the recent finite temper- 
ature lattice gauge simulations for a pure Yang-Mills SU (3) gauge theory at the known 
critical temperature. We extend this comparison using the high precision lattice data for 
two light dynamical quarks. The investigation of these three cases show some interesting 
differences arising from the strong interaction alone and in the presence of quarks. In 
this context we discuss some newer simulations for heavier quarks and other properties 
related to gluon condensation. 
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1 Introduction 



The present discussions of a possible experimental verification from heavy ion collisions 
for a new phase of matter, the Quark- Gluon Plasma (QGP), bring an added interest in 
physics at ultrahigh temperatures. This possibility relates very closely the ideas of phase 
transitions with the condensation of gluons and quarks at finite temperature. 

In this work we discuss the consequences of the gluon condensate at finite temperature 
for several different special cases. First we indicate the results of a simple calculation for 
a gluon gas with no interactions having a specified gluon vacuum expectation value as an 
ideal gluon condensate where all the interactions are in the ground state. Then we use the 
recent high precision lattice results [1,2] for the equation of state in pure SU(N C ) lattice 
gauge theory. The essential relationship for the computation of the gluon condensate at 
finite temperature [3] from the lattice data is the conformal or trace anomaly which here 
is due to the scale variance of quantum chromodynamics (QCD). It relates the trace of 
the energy momentum tensor to the square of the gluon field strengths G 2 through the 
renormalization group beta function. Here we shall discuss the approach investigated 
in [4] , for which the consequences of the finite temperature lattice data for pure SU (N c ) 
gauge theory for the gluon condensate [5] have been presented. Now, however, we are 
able to make a comparison of these previous results with some newer finite temperature 
simulations for light dynamical quarks [6], which could provide some new insights into 
the thermodynamical structure of the QGP. 

We now briefly discuss the condensation of an ideal gas of gluons, which we simply 
write as G 2 (T) for a finite temperature T. The vacuum expectation value of the pure 
gluonic system we denote by Gq, which has a known value at zero temperature [7]. From 
dimensional considerations of the structure of the gluon condensate we can write down 
the ideal gluon condensate 



where To is the temperature at which the condensation ceases. This equation gives a simple 
relation between the condensate at a finite temperature T and that at zero temperature. 
We shall soon relate G% and T to the corresponding quantities in QCD. Although this form 
for G 2 (T) as an ideal gluon gas seems presently greatly oversimplified, earlier analyses 
in finite temperature QCD had assumed [8] this behavior as an approximation to the 
decondensation of gluons. The extent to which this assumption is valid we are now able 
to state. Furthermore, we will make no attempts at the inclusion of the thermal properties 
of the hadrons in the thermodynamics of the gluon condensate. 

We start with a brief discussion of the conformal anomaly at finite temperature as it 
relates to the gluon condensate. Thereafter we look at the lattice data for the cases of 
pure lattice gauge theory and that with dynamical quarks. Finally we conclude with a 
discussion of some newer work related to these considerations, which further accent the 
fundamental differences between these three cases. 




(1.1) 
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2 The Conformal Anomaly at finite Temperature 



The relationship between the trace of the energy momentum tensor and the gluon con- 
densate has been studied for finite temperature by Leutwyler [3] for the investigation of 
problems in deconfinement and chiral symmetry. He starts with a detailed discussion of 
the trace anomaly based on the interaction between Goldstone bosons in chiral perturba- 
tion theory. Central to his discussion is the role of the energy momentum tensor, whose 
trace is directly related to the square of the gluon field strength. It is important for us 
to note that the temperature dependence of the energy momentum tensor T^ U (T) can 
be separated into the zero temperature or confined part, Tq U , and the finite temperature 
contribution #^(T) as follows: 

T^(T) = T<f + 0""(T). (2.1) 

The zero temperature part, Tq U , has the standard problems with infinities of any ground 
state. In what follows we shall not concern ourselves with the confined part since we are 
only interested here in the thermal properties. The finite temperature part, which is zero 
at T = 0, is free of such problems. We shall see in the next section how the diagonal 
elements of 9^ U (T) are calculated in a straightforward way on the lattice. The trace 0fi(T) 
is connected to the thermodynamical contribution to the energy density e(T) and pressure 
p(T) for relativistic fields and relativistic hydrodynamics [9] 

6»(T) = e(T)-3p(T). (2.2) 

The gluon field strength tensor is denoted by G% u , where a is the color index for SU(N). 
The basic equation for the relationship between the gluon condensate and the trace of the 
energy momentum tensor at finite temperature was written down by Leutwyler [3] using 
the trace or conformal anomaly in the following form: 

(G 2 ) T = (G\ - (^)t, (2.3) 
where the gluon field strength squared summed over the colors is 

-M. 

2g* 

for which the brackets with the subscript T represent thermal expectation value and the 
zero stands for the value at zero temperature. The renormalization group beta function 
/3(g) in terms of the coupling may be written as 

Thus because of the renormalization group beta function (3(g) the quantization of the 
gluon fields leads to a finite value for the average of the field strength squared at finite 
temperature which is related through the conformal anomaly to the trace of the energy 
momentum tensor. The above thermal expectation values for the field strength squared 
and the trace of the energy momentum tensor we shall write as G 2 (T) and 0Z(T) as well 
as Gq for the vacuum expectation value. 



G 2 = -^G^G%, (2.4) 
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3 Gluon Condensation in pure SU (N c ) Lattice Gauge 
Theory 



The numerical computation of the gluon condensate at finite temperature can be carried 
out in terms of quantities which are evaluated using lattice gauge theory [10]. Following 
the usual approach in in statistical physics we start with a partitition function Z(T, V) 
for a given temperature T and spatial volume V. From this basic quantity we may define 
the free energy density as follows: 

/ = ~lnZ(T,V). (3.1) 

The volume V is determined by the lattice size N a a, where a is the lattice spacing and N a is 
the number of lattice steps in a given spatial direction. The inverse of the temperature T is 
similarly determined by N T a where N T is the number of steps in the (imaginary) temporal 
direction. Thus the numerical simulation is done in a four dimensional Euclidean space 
with given lattice sizes N% x N T , which gives the spatial volume V as (N^a) 3 and the 
inverse temperature T _1 as N T a. In an SU (N c ) invariant theory the lattice spacing a is 
a function of the bare coupling f3 defined by 2N c /g 2 , where g is the bare SU(N C ) gauge 
coupling. Thereby this function fixes both the temperature and the volume at a given /3. 
Now we define the expectation value at zero temperature, Pq, as well as the spatial and 
temporal action expectation values at finite temperature, P a and P T , respectively, as the 
space-space and space-time plaquettes in the following form: 

Pa,r = 1 - ^ReiTriU^Ulut)) (3.2) 

for the usual Wilson action [2]. For the symmetric Wilson action we define the parts 
S) as 6Po on the four dimensional Euclidean symmetrical lattice of size and St as 
3(P a + P T ) on the asymmetrical lattice N 3 x N T . Now we may proceed to compute the 
free energy density by integrating these expectation values as follows: 



H i] -- .V' [cIl' S, - S T ], (3.3) 



ji4 

h 

where the lower bound (3q relates to the constant of normalization. It is important to 
note that the free energy density if a fundamental thermodynamical quantity from which 
all the other thermodynamical variables can be gotten. 

In order to compute the temperature dependence of the gluon condensate from the 
lattice data, we need the dimensionless interaction measure A(T) [10]. For this reason w 
define the lattice beta function in terms of the lattice spacing a and the coupling g as 

da- 2 

(3(g) = -2N c a^L_. (3.4) 
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Now we are able to define the interaction measure A(T) as 

A(T) = N^(g) [So - S T ] . (3.5) 

The crucial part of the recent calculations [1,2] is the use of the full lattice beta function, 
(3(g) in obtaining the lattice spacing a, or scale of the simulation, from the coupling 
g 2 . Without this accurate information on the temperature scale in lattice units it would 
not be possible to make any claims about the behavior of the gluon condensate. The 
dimensionless interaction measure is equal to the thermal ensemble expectation value of 
(e — 3p)/T 4 . Thus by the equation (2.2) above the trace of the temperature dependent part 
of the energy momentum tensor [4,5], here denoted by 9fi(T), is equal to the expectation 
value of A(T) multiplied by the thermodynamical factor of T 4 , which may be written 
(after suppressing the brakets)as follows: 

6£(T) = A(T) x T 4 . (3.6) 

There are no other contributions to the trace for pure gauge theories on the lattice. At 
zero temperature it has been well understood [7] how in the QCD vacuum the trace of the 
energy momentum tensor relates to the gluon field strength squared, Gq. Since the scale 
breaking in QCD occurs explicitly at all orders in a loop expansion, the thermal average of 
the trace of the energy momentum tensor should not go to zero above the deconfinement 
transition. Thus a finite temperature gluon condensate G 2 (T) related to the degree of 
scale breaking at all temperatures, can be defined from the trace. We have used [5] the 
lattice simulations [1,2] in order to get the temperature dependent part of the trace and, 
thereby, the value of the condensate at finite temperature. Accordingly, when we take the 
gluon condensate Gq of the confined system as the assumed value [7] of 0.012 GeV 4 for 
both the ideal gluon gas condensate and the pure gauge theory to provide a stable ground 
state. For the ideal gluon gas condensate G 2 (T) we use (1.1) with the condensation 
temperature To equal to T c , above which the condensate has completely disappeared. For 
the pure SU (N c ) around and above T c we take the published data [1,2] for A(T), and use 
the equations (2.3) and (3.6) in order to obtain the gluon condensate G 2 (T) as shown in 
Figure 1. Here we can clearly contrast the difference between the ideal gluon gas where 
the condensate has no effect at temperatures above T c and the pure SU(3) gauge theory in 
both the confined and the deconfined phases. Below T c there is a quantitative difference 
between these two cases. Above T c the ideal gas condensate is completely gone. However, 
the strongly interacting gluons for T > T c arising from the creation of gluons at these 
high temperatures continue the decondensation process into the deconfined phase. Thus 
we see clearly in Figure 1 that the pure gluon condensate never approaches an ideal gas 
for all temperatures above T c reached by simulation [1,2,4,5]. 



4 Gluon Condensation with Light Dynamical Quarks 

Now we want to compute the changes in the gluon condensate due to the presence of 
dynamical quarks with relatively small masses. There have been several extensive sim- 
ulations in recent years of the thermodynamical quantities in full QCD with two flavors 
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Figure 1: The two plots show the gluon condensate for the SU(3) lattice gauge theory 
(solid curve) as compared with that of the ideal gas (broken curve). The value for T c is 
taken as 0.264 GeV in both cases, above which the ideal gas condensate is zero. 



of staggered quarks [6,11] as well as with four flavors [12]. Although these computations 
are still not yet as accurate as those in pure gauge theory, there has been considerable 
progress in the actual precision even on not so very large lattices like 12 3 x 6. In the 
case of the more recent two flavor simulations [6] the smaller quark mass shown in Figure 
2 with the open circles is expressed in lattice units am q with the value 0.0125, which 
is actually quite small (about 0.01 GeV). The heavier quark mass of 0.025 is shown in 
Figure 2 as the filled squares. The smaller mass value in terms of the lattice size a has 
the effect of raising the temperature in the simulations so that the lighter quarks are at 
higher temperature for the same coupling. In both these cases we can see that the general 
tendency of the temperature dependence of the gluon condensation just above the onset 
of the decondensation shows a sharp drop which becomes slower at higher temperatures. 
Because of the scarcity of data points above 0.2GeV we cannot extrapolate any clear 
properties above this temperature. 

The presence of massive quarks changes the thermal properties of the gluon condensate 
due to the renormalization of the masses as expressed by the lattice beta function [12]. 
Thus the trace of the energy-momentum tensor has a direct mass dependence in the 
trace anomaly [13] for the quark condensate with the renormalized (physical) quark mass 
m q . Thereby the expectation value for the trace of the energy momentum tensor (O^) 
is calculated from the sum of the expectation values of quark and gluon condensates as 
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Figure 2: The lines show the gluon condensates for SU(3) (solid) and the ideal gluon gas 
(broken) in comparison with that of the light dynamical quarks denoted by the open circles 
and the heavier ones with filled squares. The error bars are included when significant. 



m q (ijj q ijj q ) + (G 2 ), where ip q and ?Jj q represent the quark and antiquark fields respectively. 
Here for the sake of simplicity we choose two light quarks of the same mass. We are now 
able to write down an equation for the temperature dependence of the gluon condensate [4] 
including the effects of the light quarks in the trace anomaly in the following form: 

(G 2 ) T = (G 2 ) + m q (^ q ) - mSM T - (0£) T . (4.1) 

As we have discussed previously, we identify the expectation value at finite temperatures 
for the gluon condensate with G 2 (T) in Figure 2. The simulation procedure is, in general, 
similar to the pure gauge theory with the additional effects of the quark contributions. 
The computation [6] of the interaction measure A m (T) where the subscript m denotes 
its dependence upon the quark mass m q . As before, we can calculate the temperature 
dependence of the trace of the energy momentum tensor (omitting the brackets) 

9£(T) = A m (T) x T 4 . (4.2) 

This high precision data for A m (T) is used for the compuatation of the two sets of data 
points in Figure 2 mentioned above [6]. It is possible to see that in equation (4.1) at very 
low temperatures the additional contribution to the temperature dependence from the 
quark condensate with small quark masses is rather insignificant as would be expected 
from chiral perturbation theory [3]. However, in the temperature range where the chi- 
ral symmetry is being restored there is an additional effect from the term {ip q ip q ) T . Well 
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above T c after the chiral symmetry has been completely restored the only remaining effect 
of the quark condensate is that of the vacuum. This term would contribute negatively 
to the gluon condensate. Thus we expect [5] that for the light quarks the temperature 
dependence can only contribute at all around T c . Nevertheless, because of the small quark 
mass the total effect of the light quark condensate is very small for the above mentioned 
simulated [6] values (maximally less than one procent). For the case of pure SU(N) we 
know that the values just below T c of A(T) x T 4 are very small - that is, about the same 
size as m q (i() q ip q )o. Thus, to the extent that the chiral symmetry has not been completely 
restored, its effect on G 2 (T) will be present but small just below T c . 

Finally as a conclusion to this discussion of the gluon condensate with two light dy- 
namical quarks we shall mention a few related points. While for the simulations of the 
pure SU (N c ) gauge theory we could depend on considerable precision in the determina- 
tion of T c and A(T) as well as the numerous other thermodynamical functions, this is 
still not the case for the theory with the massive dynamical quarks. The statistics for 
the measurements are generally smaller. The determination of the temperature scale is 
thereby hindered so that it is harder to clearly specify a given quantity in terms of T. 
Thus, in general, we may state that the accuracy as well as the number of data points 
for the case of the dynamical quarks are much less as compared with the computations 
of the pure lattice gauge theories. However, there is a point that arises from the effect 
that the transition temperatures with massive quarks are generally considerably lower, so 
that A(T) x T 4 is much smaller [5]. Nevertheless, there could be an indication of how the 
stability of the simulations with light dynamical quarks keeps G 2 (T) in positive values 
for T > T c so that amount of decondensation is less than the value 0.012GeV A [7] in 
the temperature range of the numerical data [6] . A model of the thermodynamics of the 
chiral restoration transition [14] has previously considered separately the condensates for 
the electric and magnetic field strengths squared. A general qualitative agreement of the 
local subtraction of the magnetic and electric condensates can be seen with our results 
in Figure 2 even though the lattice measurements which they used were not computed 
using a non-perturbative method [15], nor was the temperature scale obtained from the 
full non-perturbative beta-function [1,2]. 



5 Discussions and Conclusions 



In this work we have compared the effects of gluon condensation at finite temperature for 
three different cases- that of noninteracting gluons, the pure SU (3) thermodynamics and 
QCD with two light quarks. We have seen that the ideal gluon gas has the typical struc- 
ture of nonrelativistic Bose-Einstein condensation, whereby the condensate completely 
vanishes at a given temperature To and is exactly zero at all higher temperatures. For 
the condensate at zero temperature Gq we use the known vacuum expectation value [7]. 
This behavior we have contrasted with that for the decondensating of gluons found from 
the pure 577(3) lattice data, for which the decondensation does not stop at the critical 
temperature but continues as the temperature increases since there is no bound on the 
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number of thermal gluons [5]. The presence of two light dynamical quarks further changes 
the process of decondensation by providing other thermal mechanisms through which the 
gluon condensate begins to disappear at a much lower temperature. Although we have 
seen that this effect is slightly dependent on the value of the mass, it is strongly dependent 
on the fermionic structure relating to the chiral symmetry restoration at temperatures 
lower than the critical temperature T c of the pure SU (N c ) . Thus we are seeing the ther- 
mal effects of two different anomalies [16] at work on the same system. 

We now mention some other work which is relevant to our discussion. The very recent 
calculations performed for heavier quarks [17] provide some new aspects to the above 
analysis. At the present the pressure has been computed as a function of the temperature 
up to about 4T C for flavor numbers 2,2 + 1 and 3, where 2 + 1 means two lighter and one 
heavier quark flavors. The simulations present for other quantities considerable difficulty 
with the rriq^qip^T term when m q is larger as is the case for the s-quark. Furthermore, 
there has been some recent work on the temperature dependence of the gluon condensate 
derived from the QCD sum rules [18]. One can see from their results pronounced differ- 
ences in the thermal behavior between this method and our calculation using the lattice 
simulations [4,5]. In fact in some cases they find that the ratio G 2 (T)/Gq is larger than 
unity and growing for increasing temperatures. As a last discussion we would like to com- 
ment on some results based on the postulation of a simple form for the entropy [19]. With 
the assumption of the proper parameters relating to a bag type of model they calculate 
a form of A(T) which can lead to a similar shape of the G 2 (T) for the pure gauge theory 
shown in Figure 1. However, all of the parameters need to be determined from the lattice 
simulation or the phenomenological models independently. 
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